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Abstract 



At existing and planned neutrino factories (high energy and high intensity 



^S| ' neutrino beam facihties) precision studies of QCD in neutrino-nucleon inter- 



actions are a reahstic opportunity. We investigate charmonium production in 

^j , fixed target neutrino experiments. We find that J /tp production in neutrino- 

^^ \ nucleon colhsion is dominated by the color octet "^Si NRQCD matrix element 

,-S^ I in a neutral current process, which is not accessible in photo or leptoproduc- 

ii' tion. Neutrino experiments at a future Muon Collider will acquire sufficient 

^~^' event rate to accurately measure color octet matrix element contributions. 

,^ , The currently running high energy neutrino experiments NOMAD and NuTeV 

could also observe several such events. 
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I. INTRODUCTION 

One of the most promising future high energy facihties is the recently proposed Muon 
Colhder. In order to facihtate a decision on whether and how it should be built, all the 
various uses it can be put to should be assessed. The highly collimated and intense neutrino 
beams unavoidably generated by muon decay provide a unique opportunity for precision 
studies of QCD and electroweak physics. An excellent example of such investigations is 
related to the ongoing issue of the validity of the Nonrelativistic QCD (NRQCD) expansion 
for charmonium states and the extraction of the so-called color octet matrix elements. The 
fact that heavy quarkonium represents a non-relativistic quantum-mechanical system sig- 
nificantly simplifies its theoretical studies. In particular, the presence of several important 
scales in the quarkonium system, M , Mv and Mv"^ (~ Aqcd) where f is a small parameter 
(relative velocity of quarks in the quarkonium state) allows separation of physical effects 
occurring at different scales [jl[. 

A large excess of prompt J/ip's and ip'^^ ^^ the Tevatron over to the predictions of the 
color singlet model, i.e. the model which postulates that only quarks in a relative color singlet 
state can evolve into a charmonium, sparked both experimental and theoretical interest and 
resulted in the realization of the importance of contributions generated by the operators 
involving quark states in a relative color octet configuration. The emerging effective theory 
(NRQCD), systematically describing these processes, factorizes the charmonium production 
cross section in the form 

aiA + B^H + X) = j:-^, (0|O„^|0), (1) 

n '""c 

where Fn are short- distance coefficients containing the perturbatively calculable hard physics 
of the production of a [cc] system at almost zero relative velocity (like jg -^ ccg, qq -^ ccg, 
etc.), expressed as a series in as{mc). Here, the index n incorporates a spectral decomposition 
of the quarkonium state in terms of the quantum numbers 2S'+i^icoorj ^^ ^j-^g |^^^j gyg^^gjji^ as 
well as the number of additional derivatives acting on the heavy quark fields. The essence of 
NRQCD is to organize the above expansion in powers of the heavy quark velocity v within the 
hadron, and it can be further generalized to include other heavy quarkonium-like systems, 
such as heavy hybrids [^0. Eq. (|I]) puts all nonperturbative long-distance information 
into the NRQCD matrix elements, which describe the evolution of the [cc] system into a 
charmonium plus soft hadrons - a process that cannot be calculated at present from first 
principles. 

Several attempts have been made to determine these NRQCD matrix elements from 
various experiments. The processes involved are sensitive to various linear combinations of 
NRQCD matrix elements. The problem is aggravated by the usually very large theoretical 
uncertainties involved in these calculations (on the order of 50 — 100%), due to higher 
twist effects, uncalculated and/or incalculable higher order perturbative and nonperturbative 
contributions. In this situation any independent determination of these quantities should 
be welcome. 

A major advantage of using the neutrino beam is that, at leading order in a^, the spin 
structure of the uZ coupling selects a certain combination of octet operators. The largest 
contribution is from the one with the quantum numbers ^S"} . Of course, order of magnitude 



measurements of the size of the matrix elements of this operator have already been performed 
for the J/ip and ip' , as well as for the Xcj states. The estimates of these matrix elements 
mostly come from Tevatron fits to hadroproduction cross sections for the J/ip and Xcj and 
yield, with large theoretical errors 0], 

(0|C8/'^(=^^i)|0) ~ 0.01 GeV^ and 

(0|C>f°(35i)|0) ~ 0.01 GeVl (2) 

These values are consistent, within a ~ 50% accuracy level, with the value found from Z 
decay at LEP |Q (the latter does not separate cascade and direct production, so the value 
of 0.019 GeV^^ is understandably larger than the one in Eq. (H)). 

There are, however, large discrepancies between the Tevatron fits and the values of 
Xcj matrix elements obtained from B decays 0, and between various determinations of 
{Q\0^ (^ Si)\Q) from the Tevatron fits. Clearly, new results from HERA leptoproduction 
experiments would not clarify the situation as at leading order the process ^9 ~^ [cc\8^Si) 
is forbidden by parity conservation of strong interactions. In this situation other determi- 
nations are welcome and desired. 

The present paper is an exploratory investigation of the main features of inclusive char- 
monium production in vN collisions. This process parallels J /ip leptoproduction, in which 
case Fleming and Mehen ||^ found that the 0{a1) contribution to the total eN —>■ e J/ijj X 
cross section is small compared to the color octet 0{as) contribution. A set of cuts, requir- 
ing an energetic gluon jet well separated from the J/ip, enhances the 0{al) contributions, 
but then the color singlet contribution will dominate. These cuts, however, leave behind 
only a small part of the total cross section. We don't expect that either the difference in the 
spin structure or the (mw/mz) ~ 0.6 suppression of neutral current (NC) versus charged 
current (CC) events can change this picture, so that we feel justified to calculate only the 
0{as) contributions. We will find, however, that while the leptoproduction of J/ip is not 
sensitive to the ^S*} matrix element, and measures one combination of ^5*0 and ^PJ , 
measuring the Q^ distribution in our process allows for a determination of both the ^Sl 
and the ^PJ matrix elements. The difference is due to a difference in the spin structure of 
the Z and photon couplings. 

The relative size of the ^S*} and the ^Pj contributions to the differential cross sections 
turns out to change drastically around Q^ ~ 10 GeV^. This fact and the less steep decrease of 
the differential cross section (compared to leptoproduction) together allow an easy separation 
of the two contributing matrix elements. 

We will find that the hard process at small Q^ favors the ^S*} contribution by a significant 
factor. In an experiment which cannot easily separate direct J/ip production from a Xcj 
cascading down to a J ftp (such is the situation at LEP, for example), these cascades enhance 
the observed cross section by a significant amount. In particular, given the estimates of 
Eq. (I), we find a{uN -^ Xc2 + X) ^5a{uN -^ J/ip + X) ! Noting that BR{xc2 ^ J/^7) = 
13.5% we conclude that J /ip production via this cascade mechanism is of the same order of 
magnitude as via the direct route. It is however easy to include the effect of these unresolved 
cascades into the formalism by simply replacing the actual (OjOg (^S'i)|0) matrix element 
with a (OjOg (^5*1)10) which takes into account this cascade factor ^. The rate of XcJ 



production can then easily be related to our ^5} contribution to the J /il) production rate, 
because in that case there are no other competing matrix elements. 

The large number of events expected at the muon collider will certainly allow to study 
the entire spectrum of charmonium states. Due to the fact that the dominant ^Sq matrix 
element of the he and the dominant ^S*} of the Xcj are related by heavy quark symmetry, 
our results can be trivially translated to a prediction of these rates. The only relevant matrix 
element above is measurable in any other cascade J /ip production process. 

In the present study we also discuss what we can learn about these matrix elements at 
existing neutrino facilities. As we show below, the mere event of the detection of charmo- 
nium states in current experiments (NOMAD or NuTeV) would imply the presence of the 
color-octet structures predicted by the NRQCD factorization formalism. The main question 
to be addressed here is the smallness of the event number. Since the detectors used in cur- 
rent neutrino experiments are optimized for the observation of neutrino-related phenomena 
(such as neutrino oscillations) and not for charmonium detection, we shall concentrate on 
J/ip production, which has the cleanest experimental signature of all charmonium states. 
However cascade processes should also be included, helping to increase the event number. 
We will find that both at NOMAD and NuTeV, the event rate is on the verge of observability. 

The rest of the paper is organized as follows. In Sec. |I| we explain the velocity scaling of 
the NRQCD matrix elements that we use and provide analytical formulas for the structure 
functions F„ . We discuss the numerical results for the Muon Collider, NOMAD and NuTeV 



in detail in Sec. |T| and offer concluding remarks in the Conclusions. 



II. THE CALCULATION 

At leading order, 0{as)-, only neutral current gZ interactions contribute to charmonium 
production and Fig. [l| shows the relevant Feynman diagrams. At 0{a1), in addition to 
the color singlet contribution, one may expect some enhancement from charged current 
interactions due to the difference in the couplings and the propagators, but we do not expect 
that this could drastically change the cross section to invalidate our order-of-magnitude 
estimates. 
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FIG. 1. The leading 0{as) graphs. Note the absence of any W exchange graphs 



Before proceeding with the calculation, let us briefly review the velocity counting rules 
of the charmonium production matrix elements which will be used below. The charmonium 
production matrix elements are defined in NRQCD as 



X rrij 



(3) 



where the operators /C^'^ are bihnear in the heavy quark fields, {mj is the charmonium spin 
z-component.) Each matrix element is proportional to a power of the relative velocity v 
of the heavy quarks. Each spatial derivative acting on a heavy quark field introduces one 
power of V. The nonperturbative evolution of the cc system proceeds through multipole 
radiation of soft gluons, which introduce additional powers of the velocity. For the J/ip and 
ip', the leading matrix element is at m^v^, 



\O(^^eS,)\0) = (0|xW|J/^ + X)(J/V^ + X|^Vx|0), 



(4) 



which clearly selects the dominant Fock state. We observe, however, that this matrix element 
corresponds to a cc system in a color singlet state, which can be produced only at subleading 
order, (9(q;^). Therefore, its contribution will be smaller than that of the velocity-subleading 
color octet configuration, such as 






(5) 






2 

i 



-D^'a^^ 



T>|J/^ + X)(J/^ + X|^t 



D 



T>|J/V' + X)(J/V^ + X|V^t 



D 



-D^'a''^ 



T'^XIO), 



T"x|0), 



where {...} represents scalar, vector or traceless symmetric tensor contraction of indices. 
Five of these configurations, namely those with the quantum numbers ^5*0 , ^S"} , ^Pj 
for J = 0, 1, 2, scale as f^. The ^Pq matrix element is negligible, of the order 0{v^^), and 
will be neglected hereafter. 

In the case of the p-wave states Xcj we encounter a similar situation. The leading matrix 
elements are again suppressed by one factor of ttg, and are color singlets at m^v^. They 
project out the dominant p-wave cc combination 
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The leading matrix elements in the velocity expansion are ^Pj ~ -m^f^ and ^S*} 
but the ^S'j configuration is produced at leading order in Ug'- 






xt^T^a 



(O|0|-(^5i)|O) = (0|xVT>|xcj + X)(xcj + X|V^VT«x|0). 



t^T^«^ 



(7) 



The velocity counting rules in NRQCD require, for the most interesting case of the J/ip, 
to calculate the contributions of the 0{v'^) matrix elements, ^5*0 , ^S*} , ^PJ (heavy quark 



symmetry requires ^Pj 



(2 J + 1) ^P, 



3p(8)^ 



We note that the same structure functions can 
be used to evaluate all contributions to Xc, for which ^S"} is leading, of order 0{v^), and to 
he, for which ^Sq is leading, also of order 0{v^). We do not calculate the hard-to-identify 
r]c, where the uncalculated ^Pq would also contribute at leading order (this matrix element 
is completely unimportant for the J/ip, as it is suppressed as 0{v^^) in that case). We also 
note that parity and J conservation forbids all interference terms in Eq. (|]). 

The calculation is quite straightforward and we quote the result in the form of the 
following structure functions hn (Q^): 



da (s, Q^ 



2 2 



3sin^2^,4^ (Q2 + mz^y 



E 



(oiaio) 



m'i 



Q2+4m2 



dx fg/N {x, Q^) hn (y, Q' 



(8) 



where s is the total invariant mass of the uN system, x is the momentum fraction of the 
incoming gluon and —Q^ is the momentum-squared transferred from the leptonic system; 
y = fi_tl!2_. Because we are doing a tree level calculation, a^ does not run; in the actual 
calculation we will choose the educated guess of evaluating a^ =^ ag {Q^ + 4ml). 
Now, our calculation gives 



h i_5(8) 



h^„{s) 



'35j» 



h 3p(8) 
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3p(B) 



h-j,p(ii) 



y,Q' 



i9vf X 6 ^9,"^^^, {y' -2y + 2) 



(g2 + 4miy 



y.Q') ={9vf x2Q' 



(Q^ + I2mlf 
(g2 + 4m2)^ 



(y' -2y + 2) 



(9) 



y,Q') ={9lrf x4Q^ 



2^ _..c)2 ^^^,g^-2y + 2) + lQ{l-y)ml 

2\ , c ^2 4 ^2 (|/2 -2y + 2)Q4 + 48(1 - y)Q^ml + 96(y2 - 2y + 2)m, 



(Q2 + 4m2)^ 



y.Q') ={9vr x-Q 



(Q2 + 4mlY 

where 9% = 2 ^^^ 9v ~ 2 [^ ~ 3 ^^^^ ©vk) are the vector and axial couplings of the c-quark. 
We have checked that an appropriately modified version of these structure functions correctly 
reproduces the eN —* e J/ip X cross section in Ref . . We immediately observe that the 
coupling constants favor the '^S'} contribution, which is due to the large axial coupling (a 
similar contribution is, of course, absent in the case of J /if) lepto and photoproduction). 
Indeed our numerical estimates will show that this matrix element dominates the total cross 
section, and also the differential cross section unless Q"^ ^ ml. At large Q^, the relative Q^ 
enhancement of the P-wave structure functions makes them dominant. 

These structure functions should eventually be incorporated in the specific Monte Carlo 
generators built for each particular detector. Our numerical estimates presented in the 
following chapters ignore the fine details of the experiments and should be understood as a 
preliminary feasibility study. 
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III. NUMERICAL RESULTS 

Due to its clean experimental signature, one of the most important charmonium pro- 
duction processes is J/iIj production. We have calculated the total cross section and the 
differential Q^ distribution of J/ip^s at various experiments, both currently running and 
proposed. For lack of knowledge of the size of the NRQCD matrix elements we set them 
equal to a reasonable "baseline size" value, defined as 

{0\Os ('So) \0),aseHne = {0\Os (' S^) \0),aseHne = ^"'^'i!^'^,^^""'^ = lO"' GcV (10) 

These values satisfy the restrictions imposed by heavy quark symmetry and are compatible 
with the existing extractions @J|,5|. The results below were found using MRST parton 
distribution functions (PDF) [^, substituting rric = 1.35 GeV for the charm quark mass and 
imposing a cut Q^ > {1.2GeVy in the calculation of the total cross section. 

The energy dependence of the total cross section, taken at various representative incident 
neutrino energies, is shown in Table |. We find that the total cross section is very sensitive 
to the neutrino beam energy E,,. The origin of this strong dependence can be traced to the 
fact that the main contribution to the integral in Eq. (|]) comes from the smallest available 

X ^ — -, so that larger s probes a smaller region of x where the gluon PDF's are 

enhanced. It becomes obvious that in order to have a reasonably large event number, the 
neutrino beams with the highest possible energy should be employed. 



Eu[GeV] 


7.5 


25 


120 


450 


a[nb] 


7.8 X 10-^3 


6.9 X 10"^° 


1.3 X 10"*^ 


5.5 X 10"^ 



TABLE L Total cross sections for the J/^p production in uN -^ J/ijjX for various representa- 
tive incident neutrino energies. The values are taken to correspond to the minimal, average and 
maximum neutrino beam energy at NOMAD and maximum available neutrino beam energy at 
NuTeV respectively. 

We have calculated the total cross section for the two presently running high energy 
neutrino experiments, NOMAD at CERN and NuTeV at Fermilab, as well as for the much- 
discussed high energy muon collider. Table |I| shows the results for each term in the spectral 
decomposition. The values support our previous expectations that the ^S"} contribution 
dominates the total cross section. Note that these numbers refer to direct J/ip production; 
a cascade mechanism involving Xcj decay into J/ip can be easily calculated from our ^S"} 
contribution and its ratio to direct ^5*} production is universal and can be taken from other 
experiments. 

We should emphasize that due to the expected small number of events the question of 
background suppression becomes very important. In particular, in the case of the currently 
running experiments, NOMAD and NuTeV, the leptonic decay channel of J/ip becomes 
virtually the only possibility of detecting the produced J/ip^s. This makes electromagnetic 
lepton pair production a very important source of background [O] . 
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decomp. 


a [10-9 nb] 


10^^ 


a [10-9 nb] 


10^^ 


a [10-9 nb] 


106^ 


a [10-9 nb] 


106 ^ 


'So{8) 


0.086 


0.0131 


1.59 


1.94 


2.01 


3.89 


1.54 


2.30 


'S^i8) 


0.552 


0.0842 


9.14 


11.2 


11.5 


22.3 


8.97 


13.4 


'Po{8) 


0.129 


0.0195 


1.92 


2.34 


2.40 


4.65 


1.89 


2.82 


=^Pi(8) 


0.082 


0.0124 


2.06 


2.52 


2.64 


5.10 


1.96 


2.94 


=^P2(8) 


0.166 


0.0253 


2.71 


3.30 


3.42 


6.61 


2.65 


3.97 


SUM 


1.01 


0.154 


17.4 


21.3 


22.0 


42.5 


17.0 


25.4 


Total 


a*"* = 6.56p6 


cr*°* = 0.82 pb 


0-*°* = 0.52 p6 


a*°* = 0.67 pb 



TABLE II. The uN -^ J /ipX cross section contribution from each matrix element (set equal to 
its baseline value) in the four discussed experiments, folding in the neutrino energy distributions. 
Here, 0"*°* refers to the total deep inelastic cross section, 0"*°* = ace + c^NC- 



Another important issue is diffractive J /if) production. It is difficult to separate, both 
theoretically and experimentally, diffractive and color octet contributions. This is due to 
the fact that both of these processes contribute at z = {P^^ ■ P^)/{Pn ■ g) ~ 1. In principle, 
the absence of a rapidity gap in the color-octet NRQCD process may help . In addition, a 
perturbative QCD calculation of the diffractive leptoproduction ||ll[ suggests that this pro- 
cess falls off at high Q^ as ~ l/Q^. We expect this behavior to hold for the case of diffractive 
"neutrinoproduction" as well. As indicated in our calculation, the color octet contribution 
falls off only as ~ 1/Q'^ which makes the diffractive process negligible at sufficiently high Q^. 
Figs. 1^ and ^ show the Q^ distribution of the events. As expected, these decrease quickly 
with increasing Q^. The decrease, however, is not so drastic as not to let us use higher Q^ 
cut up to OllOGeV"^) which may also allow to reduce the nonperturbative background. 

Now we turn to the discussion of some of the presently running and planned neutrino 
experiments. 



A. The muon collider 



A proposal to build a high energy fj^^fi^ storage ring has recently drawn much attention. 
Although seemingly a sideline, luminous neutrino beams are unavoidably generated by the 
decaying muons, and provide an ideal environment for charmonium generation. Such ma- 
chines comprise both necessary ingredients, high luminosity and high neutrino beam energy, 
that are needed for detailed studies of charmonium production. 

The proposed First Muon Collider [Q would use 2 x 250 GeV muons, boosted in linear 



accelerators. These muons and antimuons decay into the i/^'s (or z7^'s) in both final booster 
(called RLA3) and in the accelerator ring. A straight section of the ring was proposed in 
order to achieve a highly collimated neutrino beam. The resulting neutrino spectra that 
have been calculated by Harris and McFarland |]13| peak in the 150 — 200 GeV region. With 
a representative set of parameters O{10^) DIS events/(g/cm^ year) are expected. This 
high number of events allows a measurement of the partial distributions without the need 
for extremely costly dedicated detectors. With the cross sections given in Table |1|, this 
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FIG. 2. The Q distribution E„[ „"^2 ] of the J/V' events at the FMC for basehne-sized matrix 
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translates to 20 — 50 J/ip events/ (g/cm^ year). 

We will now sample the various suggestions for the neutrino detectors to get a rough 
idea of the yield and capabilities. In all cases we look at a detector diameter of 40 cm, in 
which almost all of the neutrino beam fits. 

As an obvious minimal suggestion, we might look at a 30 cm thick, table-top sized de- 
tector [T^, containing approximately 50 kg of water-density target. We estimate a yield of 
1600 J /iIj events/ year, providing just enough statistics for a crude estimate of the octet 
matrix elements. By contrast, a light target considered e.g. in |T^ (aim thick liquid hy- 
drogen detector) would only produce approximately 100 J /ip events/year, not much better 
than presently running experiments. 

The best answer may be the "general purpose detector" suggested by B. King |15|. This 
detector consists of a one meter long stack of silicon CCD tracking planes, providing ~ 50 
g/cm^ density. In such a detector an estimated 3, 000 J/ip events/year would occur and 
the high efficiency and precise reconstruction capabilities should allow for disentangling the 



contribution from the color octet ^5*} and '^Pf^ matrix elements, through the measurement 



of the differential Q^ distributions. Note that an increase in the length of the straight section 
of the storage ring provides a relatively cheap way of increasing the luminosity of the neutrino 
beam, resulting in better statistics. Ref. Il^ actually considers a 200 m long straight section. 



gaining a factor of 20 over our estimates. 



Another option is a conventional fixed target type heavy detector ||T3[. A representative 



example of a 2 ton iron calorimeter 1.6 m long, would see very high event rates (approxi- 
mately 70, 000 J /ip events/year before cuts and detection efficiencies are included). If back- 
ground difficulties can be overcome, this option provides an excellent opportunity to study 
all differential rate distributions with precise statistics. 

As Fig. ^ shows, the differential cross section quickly decreases with Q^. However, the 
decrease is not as fast as in the leptoproduction case. Comparing the Q^ behavior of formula 



3c-(8) 



for the S\ contribution to Eq. (3) in Ref. 0], we find 



da 



oc 



uN 



da 



X 



Q\Q^ + (3m, 



eN 



rrir 



(11) 



where /3 ~ 4 — 10 is slightly Q^-dependent. We have checked that our curves actually repro- 
duce this relationship with j3 ^ Q. The upshot is that the absence of a photon propagator 
results in a much wider tail of the Q^ distribution than in leptoproduction, allowing for 
better discrimination between S- and P-wave contributions. 
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FIG. 3. The contributions of the various matrix elements as a function of Q^ for FMC-Ring, 
with baseline-sized matrix elements. Observe that the contributing linear combination of matrix 
elements drastically changes with Q^ within a region where the rate is still sizable, so one can real- 
istically extract both the '^S\ and the ^Pj matrix elements even without requiring exceedingly 
high statistics. 



Fig. ^ shows the ratio of the contributions from these matrix elements as a function of 
Q^. At low Q^, where the contribution into the cross section is large, the ^Si dominates 
but a determination of the differential cross section in the range ~ 2 GeV < Q < 4 GeV, 
where the event number is still sizable, would allow a reliable extraction of both ^S*} and 
'^Pj matrix elements. Higher Q^'s will only measure a linear combination of these, as the 
ratio in Fig. ^ becomes Q^-independent. 

Given the above calculation, it is easy to estimate the production rate of other charmo- 
nium states, for instance r?c. Using heavy quark symmetry, it is possible to relate, at leading 

:8) 



order in v^, the leading O^v"^) matrix elements oi rjc, ^Sq , ^Sl , ^P{ to the corresponding 
J/ip matrix elements. A simple calculation shows, using the data in Table |I| with baseline- 
sized matrix elements, that the rjc yield will be one third of the J/ijj^s. This fact and the 
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absence of a clear signature will make rjc quite elusive in this experiment. A similar remark 
holds for the he, whose only contributing matrix element, the 0{v^) level ^5*0 , is related 
to the Si matrix elements of the Xcj^s- 



B. NOMAD 

The NOMAD detector at CERN uses the CERN SPS neutrino beam, mainly designed for 
detecting neutrino oscillations. The average neutrino energy is small, 24 GeV, and, even with 
the huge mass of the detector, the event number is small. One would expect approximately 
2 — 3 events per year in the main detector with its fiducial mass of 2.7 tons. The situation is 
somewhat better in the Front Calorimeter, where such a search is underway, simply because 
of the larger mass of the detector .[] With a mass of 17.7 tons, and the detector characteristics 
given in Ref. [0, the yield becomes 14 — 15 J/ip events/year, which must be multiplied by 



the decay ratios and detector efficiencies to find the number of observable events. Given the 
inaccuracies of our calculation and the crude estimate of the size of the matrix element, it 
does not seem impossible that several such events would be seen. 

C. NuTeV 

Another high energy neutrino oscillation experiment has a chance of probing, as a byprod- 
uct, our process. The larger average neutrino energy at the Fermilab experiment, compared 
to CERN, results in approximately twenty times larger cross section, compensating for the 
smaller size of the detector. The Fermilab experiment NuTeV detected, during its lifetime 



^The authors thank Kai Zuber for drawing their attention to this point. 
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of one year, about 1.3 million DIS events. Using the numbers in Table ||, this would imply 
28 J /ip events. This is again on the verge of observability. 



IV. CONCLUSIONS 

In conclusion we note that an analogous calculation can be done for the hb system. In that 
case, however, the a^ suppression of the ratio of the color singlet vs. color octet contributions 
is compensated by the smallness of the heavy quark velocity. Indeed, we expect 

{s^nglet/octet\,_^aJv\,_^^ because ^>!lL ^ 2 ^^^ (i\\ (^_^^ ^ g. 



{singlet/ octet) cc («s/f'^)cc as{mj/^) 3' \v^ J VO.l 



(12) 



(Recall that in J/ip leptoproduction the color singlet contribution was suppressed by ~ 
10%.) As a consequence, our tree level calculation is unable to accurately predict the 
dominant T production rate. However, because all numerical factors in the color octet 
contribution are smaller for the T than for the J/ijj (including the size of the phase space), 
we can at least exclude the possibility of T production at presently running experiments. 
In order to make a meaningful estimate for the muon collider, one needs to perform much 
more involved 0{a'^) calculations. 

To summarize, we have investigated charmonium production at some of the currently 
running (NOMAD and NuTeV) and proposed (First Muon Collider) neutrino factories. 
We found that the J/tp production cross section is dominated by the contribution from 
the matrix element of the color octet ^Sl operator, which is not accessible in photo and 
leptoproduction experiments. Our exploratory study shows that the neutrino beam from the 
yU+yU" collider can be successfully used to accurately extract the color octet ^S'j and ^Pj 
matrix elements of the J/ip. This measurement is complementary to leptoproduction, where 
^5*0 and ^Pj contribute. In order to achieve an accurate extraction, however, a detailed 
investigation is necessary, including the 0(0^) contributions. A significant deviation of the 
measured values of the nonperturbative charmonium matrix elements extracted from high 
energy up and pp experiments could be a sign of intrinsic charm in the proton. 
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